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Abstract 

Relaxation effects impose fundamental limitations on our ability to coherently control quantum 
mechanical phenomena. In this letter, we establish physical limits on how closely can a quantum 
mechanical system be steered to a desired target state in the presence of relaxation. In particular, 
we explicitly compute the maximum coherence or polarization that can be transferred between 
coupled nuclear spins in the presence of very general decoherence mechanisms that include cross- 
correlated relaxation. We give analytical expressions for the control laws (pulse sequences) which 
achieve these physical limits and provide supporting experimental evidence. Exploitation of cross- 
correlation effects has recently led to the development of powerful methods in NMR spectroscopy to 
study very large biomolecules in solution. We demonstrate with experiments that the optimal pulse 
sequences provide significant gains over these state of the art methods, opening new avenues for 
spectroscopy of much larger proteins. Surprisingly, in spite of very large relaxation rates, optimal 
control can transfer coherence without any loss when cross-correlated relaxation rates are tuned to 
auto-correlated relaxation rates. 
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1 Introduction 



The control of quantum ensembles has many applications, ranging from coherent spectroscopy to 
quantum information processing. In practice, the quantum system of interest is not isolated but 
interacts with its environment. This leads to the phenomenon of relaxation, which results in signal 
loss and ultimately limits the range of applications. Relaxation is a major road block standing in the 
way of practical quantum computing. Manipulating quantum systems in a manner that minimizes 
relaxation losses is a fundamental challenge of utmost practical importance. What is the ultimate 
limit on how close an ensemble of quantum systems can be steered from an initial state to a desired 
target state in the presence of relaxation? Until now there existed no theory that answers this 
question. This situation is comparable to the time before the fundamental limits of a heat engine 
were known: More than hundred years after the invention of the steam engine, the physical limits 
for the maximum amount of work a steam engine could produce was unclear, in spite of decades 
of advances in its design. " The theory of its operation is rudimentary and attempts to improve its 
performance are still made in an almost haphazard way" pj . Of course, the maximum efficiency of 
a heat engine is not given by the cleverness of the engineer who attempts to builds such a machine, 
but by the fundamental law of thermodynamics as captured in Carnot's principle. 

In this manuscript we derive fundamental limits on how close can an ensemble of nuclear spins 
be driven from its initial state to a desired target state in the presence of relaxation. In particular, 
we derive the maximum efficiency of polarization and coherence transfer between coupled nuclear 
spins. A premier example where such coherence transfer operations are important is nuclear mag- 
netic resonance (NMR) spectroscopy 2 . In structural biology, NMR spectroscopy is an important 
technique that allows to determine the structure of biological macro molecules, such as proteins, in 
aqueous solution. With increasing size of molecules or molecular complexes, the rotational tumbling 
of the molecules becomes slower and leads to increased relaxation losses. When these relaxation rates 
become comparable to the spin-spin couplings, the efficiency of coherence transfer is considerably 
reduced, leading to poor sensitivity and increased measurement times. Recent advances have made 
it possible to significantly extend the size limit of biological macro molecules amenable to study 
by liquid state NMR [3-6] . These techniques take advantage of the phenomenon of cross-correlated 
relaxation. Cross-correlated relaxation represents interference effects between two different ralax- 
ation mechanisms 8 . Until now it was not clear if further improvements can be made and what 
is the physical limit for the coherence transfer efficiency between coupled spins in the presence of 
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cross-correlated relaxation. 



In this letter, we give analytical expressions for this maximum achievable coherence transfer effi- 
ciency for two heteronuclear coupled spins under very general decoherence mechanisms that include 
cross-correlated relaxation. We describe the optimal pulse sequences that achieve this efficiency and 
experimental data that supports these results. We demonstrate that in the limit where the inter- 
ference effects become comparable to the uncorrelated relaxation rates, complete coherence transfer 
is possible without any loss. In the general case of cross-correlated relaxation, we demonstrate 
substantial improvement over previously known sequences in NMR spectroscopy. 



2 Theory 



We consider an isolated heteronuclear spin system consisting of two coupled spins 1/2, denoted / 
(e.g. ^H) and S (e.g. ^^N). To fix ideas, we first address the problem of selective population inversion 
of two energy levels (e.g. a(3 and (3P) as shown in Fig. 1. This is a central step in high-resolution 
multi-dimensional NMR spectroscopy and corresponds to the transfer of an initial density operator 
Iz, representing polarization on spin /, to the target state 2IzSz- 



For large molecules in the so-called spin diffusion limit 2 , where longitudinal relaxation rates 
are neglegible compared to transverse relaxation rates, both the initial term (7^) and final term 
{2IzSz) of the density operator are long-lived. However, the transfer between these two states 
requires the creation of coherences which in general are subject to transverse relaxation. The two 
principle transverse relaxation mechanisms are dipole-dipole (DD) relaxation and relaxation due to 
the chemical shift anisotropy (CSA) of spins / and S. The quantum mechanical equation of motion 
(Liouville-von Neumann equation) for the density operator p is given by 

p = TT J[-i2IzSz,p] +7r kDD[2IzSz, [2IzSz,p] + n k^csAi^^^ [Iz,p]] + ti" k^sAi^z, [Sz,p]] 

+ ^DI3/CSa[2^<S'z, +71- k%^icsA^izSz,{Sz,p\ (1) 

where J is the heteronuclear coupling constant. The rates fc^^^, ^csA represent auto- 

relaxation rates due to DD relaxation, CSA relaxation of spin / and CSA relaxation of spin S*, 
respectively. The rates fc^i^j/csA ^%d/csa represent cross-correlation rates of spin / and S 
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Figure 1: Selective population inversion of the energy levels a(3 and /3/?, corresponding to a transfer 
of polarization (A) to 21^3 z (B). 

caused by interference effects between DD and CSA relaxation. These relaxation rates depend on 
various physical parameters, such as the gyromagnetic ratios of the spins, the internuclear distance, 
the CSA tensors, the strength of the magnetic field and the correlation time of the molecular tum- 
bling j2j. Let the initial density operator p(0) = A and p{t) denote the density operator at time 
t. The maximum efficiency of transfer between A and target operator C is defined as the largest 
possible value of trace(C^ p{t)) for any time t B. (by convention operators A and C are normalized). 

The main result of this letter is as follows. The maximum efficiency of transfer between the 
operators Iz and 2IzSz depends only on the scalar coupling constant J and the net auto-corrclatcd 
and cross-correlated relaxation rates of spin /, given by ka ~ kuD + ^csA ^'c ^ ^dd/csA' 
respectively. This physical limit rj is given by 

v = Vi+e-c, (2) 
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where 



c 



1,2 _ 1,2 



(3) 



The derivation of the maximum efRciency rests on the basic principles of optimal control theory and 
the development of a new class of control systems (see suppl. material for details). The optimal 
transfer scheme (CROP: cross-correlated relaxation optimized pulse) has two constants of motion. 
If li{t) and l2{t) denote the two-dimensional vectors {{Ix){t), {Iy){t)) and {{2IxSz){t), {2IySz){t)), 
respectively, then throughout the transfer process the ratio of the magnitudes of the vectors I2 and 
li is maintained constant at rj. Furthermore, the angle 7* between li and I2 is constant throughout. 
The two constants of motion of the optimal transfer scheme determine the amplitude and phase 
of the rf field at each point in time and explicit expressions for the optimal pulse sequence can be 
derived (see suppl. material). 

We now consider two important limiting cases of this problem: 

(I) In the case when fc^ > and kc = (no cross-correlated relaxation) , the optimal efRciency 77 



is equal toyl-l-^ — ■^<1 (see yellow curves in Figs. 2 and 3) and the optimal angle 7* is 7r/2 



(II) In the limit where the cross-correlation coefficient kc/ka approaches 1, the optimal transfer 
efficiency r/ approaches 1 (see black curves in Figs. 2 and 3) and 7* approaches tt. Surprisingly, in 
this case using optimal control it is possible to transfer coherence without any loss in the presence 
of relaxation. In the limit of large relaxation rates ka, this relaxation-optimized transfer mechanism 
gains up to 100% compared to state of the art transfer schemes. 



The optimal transfer scheme is best illustrated by decomposing the initial operator Iz as a sum 
of the two operators IzSa = ^ + IzSz and IzSfj — ~ IzSz- The transverse components IxSa, 
lySa and IxSfj, lySfi relax with rates ka + kc and ka — kc, respectively. When kc/ka approaches 
I, the transverse operators IxS^ and lySp do not relax. The optimal control in this case reduces 
to selectively inverting IzSp to ^IzSfj by weak rf irradiation at the frequency (—J/2) of the slowly 
relaxing multiplet component. Such selective inversions have been performed in the past in the 
absence of cross-correlated relaxation ^l^j. However, since the component IzSa, which we do not 
want to invert, has a large transverse relaxation rate given by ka + kc, it is possible to carry out the 




0. 
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Figure 2: Optimal trajectories of the two mutliplet components (A) ISa and (B) ISp for ka = J 
and kc/ka = (yellow curves), kc/ka = 0.75 (red curves), kc/ka = 0.95 (blue curves), kc/ka = 0.999 
(black curve). 

selective inversion process much more rapidly. In Fig. 2, optimal trajectories of the two multiplet 
components are shown for several cross-correlation coefficients kc/ka and ka — J. 

In Fig. 4, the optimal rf amplitude and irradiation frequency of a CROP sequence is shown 
as a function of time for the case kc/ka = 0.75 and ka — J- Although the ideal sequence has a 
long duration, most of the transfer occurs in a relatively short time window, outside of which the rf 
amplitude is vanishingly small. The transfer efficiency rj is shown in Fig. 3 A for several ratios kc/ka 
as a fimction of the auto-correlated relaxation rate ka/J- For the case kc/ka — 0.75, the physical 
limit of the transfer efficiency is compared in Fig. 3 B to the transfer efficiency of conventional 
transfer schemes. 

The optimal control methods for the transfer from to 21 ^Sz in the presence of cross-correlated 
relaxation immediately extend to other routinely used transfer, such as inphase to inphase transfer 
{Ix — > Sx) |S| and single transition to single transition transfer {21 xS" — > 21" Sx) 4_. Since the 
operators Iz, Sz and 2IzSz do not decay, the optimal efficiency for the transfer Ix to Sx is achieved 
by first rotating Ix to Iz (which can be done rapidly with neglegible loss). Then 1^ is transferred 
optimally to 2IzSz with efficiency rj (Eq. 2), followed by the optimal transfer of 2IzSz to Sz, which is 
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finally rotated rapidly to Sx- The optimal transfer 21^3 z ~> Sz is analogous to the optimal transfer 
Iz 2IzSz- The efficiency 77' for this transfer is also given by Eq. (2), where the rates ka and 
are replaced by the corresponding rates k'^ ~ kon + f^csA ^^^"^ ~ ^dd/csa ^'^^ ^P^-"^ ^ ^ 
replaced by the corresponding The maximum efficiency for the transfer — > Sz is the product 
of the efficiencies of the individual steps (see Table 1). 
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Figure 3: (A) Physical limits of the transfer efficiency as a function of ka/J for k^/ka = (yellow 
curve), kc/ka = 0.75 (red curve), kc/ka = 0.95 (blue curve), kc/ka = 1 (black curve). (B) For the 
case kc/ka = 0.75, the theoretical bound of the transfer efficiency (CROP: red curve) is compared 
to the transfer efficiency of conventional transfer schemes (INEPT: blue curve, CRIPT: green curve, 
CRINEPT: purple curve). 
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Figure 4: Truncated cross-correlated relaxation optimized pulse (CROP) for kc/ka — 0.75 and 
ka = J: Radio frequency amplitude A = ^Bq/^-k (left) and irradiation frequency v (right) as a 
function of time. 

In the light of increasing use and superiority of TROSY (Transverse Relaxation-Optimized Spec- 
troscopy) methods 4 , the single transition to single transition transfer 2/3, 5" 21°^ Sx is important 
in NMR applications to structural biology. It is of both theoretical and practical interest to establish 
the physical limits for this transfer. This transfer can be achieved optimally as a sequence of the 
following steps. First the term 2IxS" is rapidly rotated to 2IzS°' = Iz + 2IzSz- In a second step, 
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Table 1: Bounds on Coherence and Polarization Transfer 



Transfer 



Physical Limits of Efficiency 




eta' = + C'^ - C 
7777' 



X 



2IzSz is transferred via CROP to Sz-, followed by the CROP transfer of to 2/^5^. This completes 
the transfer from 21 ^S" to 21°" Sz which is finally rapidly rotated to 21°' Sx- The maximum overall 
transfer efficiency is given by x/772 + 77'^ (c.f. Table 1). 

3 Experimental results 

The performance of the analytically derived CROP sequences was tested experimentally using the 
coupled two-spin system of ^^C-labeled sodium formate with a coupling of J = 193.6 Hz between 
the ^^C spin (denoted /) and the spin (denoted S). In order to control the rotational correlation 
time, sodium formate was dissolved in a mixture of 96% Dg-glycerol and 4% D2O. The viscosity of 
this solvent can be conveniently adjusted through a variation of temperature. The experiments were 
performed at a temperature of 256.5 K where ka/J ~ 1 (see Fig. 5 A) and 260 K where ka/J ~ 0.5 
(see Fig. 5 B). At a magnetic field of 17.6 T, the experimentally determined ratio of cross and 
auto correlation rate was kc/ka ~ 0.75. In the preparation phase of the experiments, the thermal 
equilibrium magnetization was dephased by applying a 90° proton pulse followed by a pulsed 
magnetic field gradient. The transfer efficiency of ^'^C polarization Iz to 2IzSz was measured for the 
novel CROP sequence, as well as for INEPT |T21, CRIPT [m and CRINEPT sequences. Finally, 
a hard 90° proton pulse was applied to transform 2IzSz to 2IzSx and the amplitude of the resulting 
proton anti-phase signal was measured. The resulting experimental transfer amplitudes arc shown in 
Fig. 5 as a function of the transfer time. CROP sequences were truncated symmetrically to acquire 
transfer amplitudes also for finite mixing times. Experimentally, the optimal transfer time of the 
CROP sequence was found to be 7.5 ms. This is a compromise between losses due to the truncation 
of the (very long) CROP sequence and losses due to the non-zero relaxation rates of the terms IzSz- 
The experimentally determined relaxation time of these terms was about 50 ms. In spite of these 
non-idealities of the model system, the CROP sequences are substantially more efficient than the 
conventional sequences. In Fig. 5 A and B, the experimental gains compared to CRINEPT are 34% 
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and 22%, respectively. We found the although the optimal pulse sequences were designed for specific 
rates ka and fcc, they were robust to variations in these parameters. 
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Figure 5: Experimental transfer amplitude of truncated CROP sequences compared to CRINEPT, 
CRIPT, and INEPT as a function of total transfer time. The experiments were performed at a 
temperature of 256.5 K with ka/K « 1.1 (A) and 260 K with ka/kc « 0.6(B). 



4 Conclusion 



Here, we derived for the first time upper achievable physical limits on the efficiency of coherence and 
polarization transfer for two coupled spins in the presence of very general decoherence mechanisms 
that include cross-correlated relaxation. In this letter, the focus was on the study of polarization and 
coherence transfer between an isolated pair of scalar coupled heteronuclear spins in the spin diffusion 
limit. For this example, new transfer schemes were found which yield substantial gains (of up to 
100%) in transfer efficiency over conventional methods. With these physical limits established, it is 
expected that significant improvement can be achieved over state of the art experiments in protein 
NMR spectroscopy. Work is in progress to incorporate practical considerations like broadbandedness 
and robustness with respect to variations of relaxation rates and experimental imperfections. The 
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methods presented here can be generahzed for finding relaxation optimized pulse sequences in larger 
spin systems as commonly encountered e.g. in backbone and side chain assignments in protein NMR 
spectroscopy. Furthermore these methods directly extend to other routinely used experiments like 
excitation of multiple quantum coherence . The most surprising aspect of the presented results is 
that in spite of large relaxation rates, it is possible to exploit the structure of relaxation and have 
decoherence-free evolution by steering the system through a decoherence-free subspace. It is also 
expected that the methods presented here will be further developed to minimize decoherence losses 
in various proposed implementations of quantum information processing. 
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5 Supplementary Material: Derivation of Optimal Control 
of Coherence Transfer under Cross-Correlated Relaxation 



The derivation of optimal control assumes that the two heteronuclear spins have well separated 
resonance frequencies, allowing for fast selective manipulation of each spin on a time-scale determined 
by the coupling J and the relaxation rates ka and kc. Therefore the Cartesian spin operator can 
be transformed to an operator of the form Ix cos /3i + Iz sin j3i by the use of strong, spin-selective 
radio frequency (rf) pulses without relaxation losses (see Fig. 6). Let ri{t) represent the magnitude 
of polarization and in-phase coherence on spin / at any given time i, i.e. r1{t) — {Ix)"^ + {ly)"^ + {Iz)"^ 
(where e.g. {Ix) ~ tracejp Ix} represents the expectation value oi Ix)- Let li{t) be the magnitude 
of in-phase coherence on spin /, i.e. l\{t) — {Ix)^ + {ly}^ and Pi = cos~^ ^ (see Fig. 6). Using rf 
fields, we can exactly control the angle (3i. Hence we can think of cos/3i as a control parameter and 
denote it by Ui (see Fig. 6). Observe that the operator I^ is invariant under the evolution equation 
(1), whereas Ix and ly evolve under the J coupling and also relax. For example, Ix evolves under 
the coupling to 2IySz with rate J and cross relaxes to ~2IxSz with rate kc- In the plane defined by 
the operators 2IxSz and 2IySz, the direction in which antiphase coherence begins to build up from 
an initial coherence Ix forms an angle 

= tan-i(^) 

with the axis 2IxSz- 



<I7> 




<2hS,> 



y</x>'+</y>' 




y<245z>'+<2/y5z>' 



Figure 6: Representation of the system variables ri, r2, their transverse components h, I2, the angles 
(3i, ^2, and of the control parameters ui = cos/?i, U2 = cos/32 in terms of the expectation values 
(Ix), {ly), (Iz), {2IxSz), {2IySz), and {2IzSz). 
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As the operators 2IxSz and 2IySz are produced, they also relax. Let hit) measure the magnitude 
of the total antiphase coherence at time t, i.e. = {21^.3^)'^ + {2IySz)'^. By use of rf pulses 

it is possible to rotate the antiphase operators 2IxSz and 2IySz to 2IzSz, which is protected from 
relaxation. Let r2 represent the total magnitude of the expectation values of these bilinear operators, 
i.e. r|(i) = I2 + {2IzSz)^ and — cos~^ ^ (see Fig. 6). We can control the angle /?2 and we define 
cos/?2 as a second control parameter U2 (see Fig. 6). 







</y> 








<2/y5z> 
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<Ix> 
<2/xSz> 

Figure 7: The figure shows the transverse planes defined by Ix and ly superimposed on the plane 
defined by 2IxSz and 2IySz such that Ix is aligned with 2IxSz- h and I2 are the vectors representing 
transverse coherences in these two planes and 7 is the angle between the vectors. 



We superimpose the transverse planes defined by Ix and ly with the plane defined by 2IxSz and 



2IySz such that Ix is aligned with 2IxSz (see Fig. 7). If 7 represents the angle between li and h, 
then by definition of 6 we have 
d 



dt 
d 

di 



hit) 
hit) 



-7r[ kahit) - v/^fTj2 cos(0 + 7) hit)] 

-7t[ kahit) - ^kl + J2 COS(0 - 7) hit)]. 



This can be rewritten as 

d_ r hit) 

Jt [ hit) 

where 



ttJ 



-i X cos(6' + 7) 

X cos(6l-7) 



^ = ka/J 



hit) 
hit) 



(4) 



and 
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For a given value of Pi and (32, we have li{t) — ri{t) cos(/3i) and hit) ~ r2{t) cos(/32)- Since ~ 
and '^^^^t'^'^ = 0, we get ri{t) — ii{t) cos(/?i) and r2{t) = /2(i) cos(/32)- Substituting for /i and I2, 
we then get 



d_ 

dt 



ri{t) 
r2{t) 



X U1U2 cos{9 - 7) 



X U1U2 cos{9 + 7) 



nit) 

r2{t) 



(5) 



Given the dynamical system in equation (O, we want to find the optimal values of U2{t), 
and 7(i), so that starting from ri(0) = 1 we achieve the largest value for r2. 



As the operator is transferred to 21 ^Sz, the ratio ^ increases from to 00. The optimal 
choice of ui and U2 must ensure that the ratio of gain dr2 in r2 to loss dri in ri for incremental time 
steps dt is maximized at each step. This ratio is 

|r'2| _ -Cu|r2 + xuiM2''i cos(6' - 7) 
k'l I C^i'"! - X"i"2 cos(6' + 7)r2 

Let = g, then the above expression can be re- written as 



ri -^g^ + x g cos{9 - 7) 



|r'i| r2 ^ - X 9 cos{e + 'j) 

This expression needs to be maximized with respect to choice of g and 7. Let these optimal values 
be rj and 7* respectively. Then 



yields 



dg 



7*) + 77Cos(6' + 7*) 



2^ 



(6) 



which yields p4 — —rf'. Now the value of 7* in equation 10 is such that it maximizes 77. Differ- 



In I 

entiating both sides of equation with respect to 7* and substituting 



0, we obtain that 



-sin(6l-7*)-77sin(6l + 7*) = 0. 

■n 



(7) 



The optimal 77 and 7* then satisfy equations 0) . The two equations can then be solved to give 



p- and optimal 7* = tan 



1-C, (8) 
(1+1)^) cot fl ' substituting the optimal control law 



where ( - y 72+^2 

'''^^ -, and integrating equation lO, we see that r2 increases from to + 1 - C, which 



U2{t) 



is then the maximum achievable transfer. 
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We see that throughout the optimal transfer, the angle 7 is maintained constant at 7*. The optimal 
control ensures that the ratio of the transverse components hit) and li{t) is always maintained 
constant at rj. These two constraints can now be used to get explicit expressions for the magnitude 
and phase of the optimal rf-field. 

Let (f) denote the phase of the rf-field relative to h { j — (f> relative to I2 ) and let A be its 
amplitude. Let dli denote the change in transverse component perpendicular to the vector li 
by application of the rf field in small time dt. Then observe dli — 2TiA{Iz)cos{(f))dt. Similarly 
dl2 = '2ttA{2IzSz) cos(7 — (l>)dt. If ^ is maintained at ry then the angle 7 does not change due to the 
evolution equation (0J. Therefore we only need to consider the change in 7, due to the rf field. If 

7 is maintained constant, then = . This gives tan(/3i) cos((/)) = tan(/32) cos(7* — 0) because 

(I ) (21 s ) 

— tan Pi and -^^ — — tan /?2 . This then implies that 

_i / tan/3i , .\ 

^ = tan — - cot 7 . 

\tan/52sm7* / 

The amplitude can be determined from the condition that ^^^jy is maintained constant. This 
implies that ^ = Substituting dh = {-£,ttJIi - x7rJcos(6' + 7)^2 + 2TrA{Iz) sm{(l)) )dt and 
dl2 = (— ^7rJ^2 — X7rJcos(6' — 7)^2 — 2ttA{2IzSz) sin(7 — 0) )dt, we get 

^_ 1 (cos(6'-7*) -7?2cos(6' + 7*))x^ 
27r (tan /3i sin + tan /32 sin(7* —(/))) 

The expressions of A and cj) are given in terms of the state of the system (angle f3i and (32)- We 
can insert these in the equations for how Pi and P2 evolve as a function of A and to get explicit 
expressions of P2, 4> and A as a function of time. This gives us the amplitude and phase of the 
optimal rf pulse as a function of time. 
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